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Flow measurements are reviewed with particular emphasis on the hydrodynamic char- 
acter of elliptic flow at RHIC. Hydrodynamic scaling compatible with the production of 
highly thermalized matter having a high degree of collective interactions and extremely 
low viscosity, is observed for a broad selection of the data. These properties suggest the 
production of a new state of strongly interacting nuclear matter at extremely high density 
and temperature where the relevant degrees of freedom are the valence quarks, ie. the 
sQGP. 

1. Prologue 

For more than twenty five years the heavy ion community has sought to create and 
study the quark gluon plasma (QGP) - a new phase of hot and dense nuclear matter 
where quarks and gluons are no longer confined to the interior of single hadrons. Such a 
state is predicted to exist at high energy densities by quantum chromodynamics (QCD) 
and is now very strongly indicated by lattice QCD calculations [HH2J- ^he fundamental 
value of this search is rooted in the fact that the QGP; (i) is the ultimate primordial form 
of QCD matter at high temperature or baryon number density, (ii) was present during 
the first few microseconds of the Big Bang, (iii) provides an example of phase transitions 
which may occur at a variety of higher temperature scales in the early universe, and (iv) 
can provide important insights on the origin of mass for matter, and how quarks are 
confined into hadrons [ 3j. 

The use of flow correlations as a probe for the nuclear equation of state (EOS) and 
a possible phase transition, was recognized quite early [ El El El Ej The connection is 
made transparent in the framework of perfect fluid dynamics where the conceptual link 
between the conservation laws (baryon number, and energy and momentum currents) 
and the fundamental properties of a fluid (its equation of state and transport coefficients) 
is straightforward. It is therefore not surprising that hydrodynamic considerations have 
played, and continue to play an important role in flow studies [Ej. Fig.[I]shows the result of 
an early fluid-dynamical calculation which illustrates the development of flow correlations 
in Ne+U collisions. Fig. El shows a schematic diagram of a more recent conjecture that 
jet interaction with a high energy density medium could lead to the generation of a Mach 
cone and hence conical flow [ 9j. 
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Figure 1. Density and temperature con- 
tours in the scattering plane of a mid-central 
Ne+U collision. The arrows indicate veloc- 
ity fields. Results are from a fluid-dynamical 
calculation [[7j. 
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Figure 2. Schematic diagram of the devel- 
opment of a Mach cone resulting from the 
strong interaction between a jet and a high 
energy density medium [ 9j. The arrows 
which point to the left indicate the direction 
of flow normal to the shock cone. 



1.1. Azimuthal distributions and flow correlations 

Azimuthal angle distributions and correlation functions play a pivotal role in the study 
of flow correlations. Experimentally, one measures the magnitude of flow correlations by 
evaluating the Fourier coefficients A n , of the anisotropy of the distribution in azimuthal 
angle difference (A0 = — 2 ) between pairs of charged hadrons [fTUl fTTl [T2] 

E \ n = (e inA *). (1) 

Here 0i ? 2 are the azimuthal angles of a pair of particles measured in the laboratory coor- 
dinate system, and the brackets denote an average over pairs of particles emitted in an 
event followed by further averaging over events. Alternatively, the Fourier coefficients v n 
can be obtained as [ El El El E| : 

v n = ( e '»(*-*K)) = (cosn(0 - $ fl )>, (2) 

where (j) is the azimuthal angle of an emitted particle (also measured in the laboratory 
coordinate system) and is an estimate of the azimuth of the reaction plane. A requisite 
correction which takes account of the dispersion of the reaction plane is easy to evaluate 

[ usi nn uni uni 

For detectors having a limited acceptance, the correlation function C2(A0), is often 
exploited to measure A n [ El El Ej : 
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where the correlated distribution N cor (A0), is obtained from pair members belonging to 
the same event. The mixed distribution N mix (A0), is made of pair members belonging 
to different events having characteristics (multiplicity, vertex, etc.) similar to those for 
N cor (A0). The correlation between each particle and the reaction plane induces corre- 
lations among all particles. Thus, it is straightforward to show that X n = v\ when both 
particles are selected from the same pt and rapidity y ) range (ie. fixed-p^ correlations) 
and \ n = v n (l)v n (2) when particle (1) and (2) are selected from a different pt and/or 
y range (ie. assorted correlations). When v n is measured for specific particle species or 
as a function of transverse momentum, centrality and/or rapidity, they are referred to as 
"differential" flow measurements. When the measurements are averaged over a sizable 
phase space they are termed "integral" measurements. The first two coefficients, V\ and 
v 2) are usually referred to as directed flow and elliptic flow, respectively. 

In addition to the consideration of detector acceptance, reliable flow measurements 
often require the suppression of so called non-flow correlations. These include small-angle 
correlations due to final state interactions and quantum statistical effects [Eh correlations 
due to resonance decays [Ej and mini-jet production [Ej- The methods of Lee Yang zeros 
[EDj an d cumulants [ IZT] have been developed to suppress such non-flow correlations. The 
cumulant method is based on a cumulant expansion of multi-particle correlations [ IZT], 
The method of LeeYang zeros follows the spirit of the Lee Yang theory of phase transitions 
[ Ej ) and extracts flow directly from the genuine correlations involving a large number of 
particles [Ej. 

Azimuthal distributions are not only important for flow measurements. They are of 
great current interest in connection with the study of jet modification [Ej, the formation 
of disoriented chiral condensates [ElEj, and the study of parity and/or time-reversal 
violation [ Ej- The combined use of flow correlations and two-particle interferometry 
measurements is also used extensively to gain detailed insight on the three-dimensional 
structure of emitting sources [ El El Ej • 

1.2. A snapshot of over twenty years of elliptic flow measurements 

Systematic flow measurements now exist for a beam energy range which spans nearly six 
orders of magnitude. They include measurements done at MSU, GANIL, BEVALAC, GSI, 
Dubna, AGS, SPS and RHIC. Results for elliptic and directed flow in Au+Au collisions, 
from several of these measurements are summarized in Figs. El and ^ respectively. For 
a broad range of beam energies (Ef, eam /A > 400 MeV), the elliptic flow results can be 
understood in terms of a delicate balance between (i) the ability of pressure developed 
early in the reaction zone, to effect a rapid transverse expansion of nuclear matter, and (ii) 
the passage time for removal of the shadowing of participant hadrons by the projectile and 
target spectators [ El El Ej- The characteristic time for the development of expansion 



perpendicular to the reaction plane is <~ R/c S) where the speed of sound c s = ydP/de, 
R is the nuclear radius, P is the pressure and e is the energy density. The passage time 
is rsj 2R/ (70^0) ? where v is the cm. spectator velocity. 

If the passage time is long compared to the expansion time, spectator nucleons serve 
to block the path of participant hadrons emitted toward the reaction plane, and nuclear 
matter is squeezed-out perpendicular to this plane giving rise to negative elliptic flow. The 
squeeze-out contribution should then reflect the ratio c s j (70^0)- This is put into evidence 
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Figure 3. Elliptic flow (v 2 ) vs beam energy Figure 4. Directed flow vs. beam energy 
for mid-central collisions Au+Au colllisions [ for mid-central Au+Au collisions. The lines 
IHlj . show transport theory predictions for differ- 

ent EOS's [ESj. 



in Fig. EH where the differential elliptic flow values v 2 (pt)i shown for beam energies of 0.4 
- 2 AGeV in Fig. are scaled by the passing time. The rather good scaling observed 
suggest that c s does not change significantly over this beam energy range. 
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Figure 5. Differential flow v 2 (pt) vs pt for Figure 6. The same data shown in Fig. El 
several beam energies. The data are from scaled by the passing time. 
Refs. [H3E3 



For very short passage times (such as those at top AGS energies and beyond), the 
inertial confinement of participant matter is significantly reduced and preferential in- 
plane emission or positive elliptic flow is favored. This is so because the geometry of the 
participant region exposes a larger surface area in the direction of the reaction plane and 
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pressure gradients in this direction are also larger. Thus, the observed trends ( ie. negative 
for beam energies < 4 AGeV and positive for beam energies > 4 AGeV) for the beam 
energy dependence of elliptic flow are well understood [ E21 ESI E2j For RHIC energies, 
strong Lorentz contraction and very short passage times lead to significant reduction of 
shadowing effects, and positive elliptic flow is driven primarily by the anisotropics in the 
transverse pressure gradients [ 381 HHj . 




Ebeam/A(GeV) p/p Q 



Figure 7. Differential flow v 2 vs ^ eam for Figure 8. Zero temperature EOS. The 
the energy range 0.4 < E^am < 12.0 AGeV. shaded region is the region of pressures con- 
Results are shown for protons for mid-central sistent with the experimental flow data. The 
collisions [ESI- curves show predictions for different EOS's [ 

ESI 



1.3. Pre-RHIC constraints for the EOS 

The elliptic and directed flow observables are sensitive to the mean field and to the 
equation of state. Consequently, comparisons of these two observables to model calcula- 
tions can provide important constraints for the EOS. Such comparisons have been carried 
out for the beam energy range 0.4 < E beam < 12.0 AGeV 

Figs. 0] and show comparisons for several incompressibility constants (K). At incident 
energies of 2-6 GeV/A, for example, the transverse flow data lie near or somewhat below 
(to the low pressure side of) the K = 210 MeV calculations, while the elliptic flow data lie 
closer to the K = 300 MeV calculations. The calculations without a mean field (cascade) 
or with a weakly repulsive mean field {K = 167 MeV) provide too little pressure to 
reproduce either flow observable at higher incident energies (and correspondingly higher 
densities). The calculations with K = 167 MeV and K=380 MeV provide lower and upper 
bounds on the pressure in the density range 2.0 < p/po < 5.0. 

If one factors in the uncertainties due to the momentum dependencies of the mean fields 
and the collision integral, a range of pressure-density relationships can be established from 
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the comparisons made in Figs. 0] and These bounds on the equation of state are shown 
for zero temperature matter by the shaded region in Fig. [HI They constitute the most 
current set of constraints on the EOS for high energy density nuclear matter created at 
high baryon densities, and can provide a rudimentary base-line for comparisons involving 
the essentially baryon- free matter created at RHIC. 

2. Elliptic flow at RHIC 

2.1. The nature of the "soft" matter formed at RHIC 
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Figure 9. Yield vs transverse energy E t . The Figure 10. (E T ) / (N c h) vs. ^s^n for 
data are obtained from Ref. [HI|. Au+Au collisions. The data are taken from 

Ref. [E3J. 



It is a well known fact that rather high energy densities are achieved in central and 
semi-central heavy ion collisions at RHIC [EJEBj. Fig. [HI shows a recent measurement of 
the transverse energy E T) distribution for central Au+Au collisions at -^/snn = 130 GeV. 
Following the commonly exploited Bjorken ansatz 

1 1 dE T 

£Bj 7r R 2 T dy 1 

if we assume a thermalization time r ^0.2 — 1 fm/c, one obtains the energy density 
£Bj = 5 — 15 GeV/fm 3 , which is 35 - 100 times the energy density £ , f° r normal nuclear 
matter. This energy density far exceeds the lattice QCD estimate (~ 1 GeV/fm 3 ) for 
creating a de-confined phase of quarks and gluons (QGP)[[IJ Ej. Fig. fTUI shows the ^/s^n 
dependence of the ratio (Et) / (A^), obtained for Au+Au collisions spanning the collision 
energy range SIS - RHIC. It shows that the (Et) rises faster than the mean multiplicity 
for charge particle production (N c h), from SIS to SPS energies. However, there is little 
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Figure 11. Comparison of the multiplicity Figure 12. Comparison of the multiplicity 
distributions for Cu+Cu and Au+Au colli- distributions for Cu+Cu and Au+Au colli- 
sions obtained for the same number of par- sions obtained for the same number of par- 



ticipants at y/SNN — 62.4 GeV. The data are ticipants at ^snn = 200.0 GeV. The data 
obtained from Ref. [ESj- are obtained from Ref. [Ex- 

change in the ratio (E T ) / (N c h) from SPS to RHIC, suggesting that the s/snn increase 
results primarily in an increase in particle production. 

There are strong hints that this particle production follows rapid thermalization. First, 
the multiplicity distribution for the same number of participants is independent of colliding 
system size, as would be expected from a system which "forgets" how it is formed. A 
beautiful demonstration of this is shown in Figs. [H] and [El with PHOBOS data for 
Au+Au and Cu+Cu collisions at two separate collision energies. Second, a comparison 
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Figure 13. Comparison of measured particle ratios for Au+Au colli- 
sions, with the predictions of a statistical model (solid lines). 
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of the measured anti-particle to particle ratios with the predictions of a statistical model 

^ e -(E+y)/T 
p e -{E-»)lT 

gives excellent agreement as shown for Au+Au collisions in Fig. El It is noteworthy that 
the the temperature T ^ 176 MeV, required for this excellent agreement compares well 
with the critical temperature T cr ^170 MeV, required for the QGP phase transition. The 
baryochemical potential /i, extracted for this temperature, increases from ~ 29 MeV at 
^/s^ = 200.0 GeV to - 40 MeV at y/s^ = 130.0 GeV. 

Given the large energy densities produced in RHIC collisions, if thermalization does 
indeed occur, then one expects the development of large pressures; 



p = p- ( d£ /d P . 



s/p 

More importantly, if thermalization is rapid in non-central collisions, then large pressure 
gradients resulting from the initial spatial anisotropy or eccentricity 



of the collision zone would lead to strong elliptic flow. That is, the observation of large 
elliptic flow for a variety of particle species would be compatible with the expectation for 
early thermalization. 

2.2. Elliptic flow results 



Au+Au 




Figure 14. PHENIX preliminary data Figure 15. PHENIX preliminary data 

for the doubly differential anisotropy for the doubly differential anisotropy 

V2{centrality ) pT) ) for mesons (dominantly V2{centrality ) pT) ) for protons and anti- 

pions). protons. 

Detailed differential and integral measurements are now available for charged hadrons 
and a variety identified particle species (p, 7r, K ) A, Q, 0, 5, d, D) at RHIC [ HZj. Figs. El 
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and [K] summarize the extracted doubly differential anisotropy v 2 (centrality , p?) , for 
mesons (dominantly pions) and baryons (protons and anti-protons). They give an ex- 
cellent overview of the the detailed evolution of v 2 as centrality and pt are varied. The 
results shown for protons (p) and anti-protons (p) give an especially good view of the 
evolution away from the well known quadratic dependence of v 2 (pt) (which is also ob- 
served in very central collisions for these data) as the collisions become more peripheral. 
Figs. [03 and El show a representative set of differential flow measurements v 2 (pt) for 
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Figure 16. Preliminary STAR data showing 
v 2 (pt) vs. pt for different particle species 
produced in Au+Au collisions at s/snn = 
200 GeV [S7j. 
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Figure 17. Preliminary PHENIX data show- 
ing v 2 (pt) vs. pt for different particle species 
produced in Au+Au collisions at s/snn = 
200 GeV [S7j. 



several different particle species. They are all characterized by rather large magnitudes 
compatible with the predictions of the hydrodynamic model which in turn 

implies the creation of a strongly interacting medium and essentially full local thermal 
equilibrium. It is especially noteworthy that the v 2 values for the (j) are comparable to 
those of other hadrons. Given its rather small re-scattering cross section, such v 2 values 
suggest that thermalization is rapid and pressure gradients develop at the partonic level. 

Further evidence for rapid thermalzation can also be found in the recent measurements 
of v 2 for non-photonic electrons [ H2j . Such measurements provide a probe for whether 
or not the charm quark flows. Fig. [H3 compares recent PHENIX measurements to the 
calculations of Rapp et al. [H21- The good agreement (for p T < 2 GeV/c) between the 
data and the calculation which includes charm flow, suggests heavy quark thermalization 
in RHIC collisions. 

As discussed earlier, the pressure or energy density gradients provide the fluid accel- 
eration necessary to develop the large observed values of v 2 . The initial energy density 
controls how much flow can develop globally, since it sets the overall time scale between 
the beginning of hydrodynamic expansion and final decoupling. On the other hand, the 
detailed development of the flow patterns are controlled by the temperature dependent 
speed of sound c 2 s = dP/de. Lattice QCD calculations indicate that c 2 s « 1/3 for T > 2T cr , 
but drops steeply by more than a factor of six near the "softest point" where T « T cr . 
Thereafter, it rises again in the hadron resonance gas phase to the value c 2 s ~ 0.15. Given 
this, an important question is whether or not there might be a trace of this "softest point" 
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Figure 18. v 2 (p T ) vs. p T for charm [H2j- The Figure 19. <s/s NN dependence of v 2 for fixed 
lines compare the results from a calculation p T selections of 0.65 and 1.75 GeV/c from 
which includes/excludes charm flow [ESI- A+A collisions with A ~ 200 [IHU]. 

in the data. 

Fig. fT9l shows the ^s^n dependence of v 2 over a broad range of bombarding energies for 
two pt selections. The data show essentially no beam energy dependence over the entire 
energy range explored at RHIC, but decreases substantially as one moves down to the 
lower SPS and AGS energies. These data reflect the predicted non-monotonic structure 
in the elliptic flow excitation function [Ell- However, they do not show the decrease (due 
to very strong softening of the EOS in the phase transition region) followed by a recovery 
in the moderately stiff hadron gas phase predicted by the model [EI]- Nonetheless, the 
apparent saturation of v 2 above s/s^n = 62.4 GeV does not exclude the role of a rather 
soft equation of state resulting from the production of a mixed phase [ E2j for the range 
= 62.4 - 200 GeV. 

3. Universal scaling and perfect fluid hydrodynamics 

A particularly important question of great current interest is whether or not the fine 
structure of azimuthal anisotropy (ie. its detailed dependence on centrality, transverse 
momentum, particle type, higher harmonics, etc) can provide valuable constraints on (i) 
the range of validity of perfect fluid hydrodynamics, (ii) the parameters of the model and 
(iii) the onset of competing mechanisms such as quark-coalescence [ 53 . One approach 
to address such questions, is the detailed investigation of flow data to test the validity 
and/or failure of several scaling "laws" predicted by perfect fluid hydrodynamics. 

An important scaling prediction of hydrodynamic theory is exemplified by the exact 
analytic hydro solutions [ 54J exploited in the Buda-Lund model [ESI ES|- The model 
gives: 



, Pt 

IqIw) 4m t 




(4) 



where 7 ,n are modified Bessel- functions, m t is a rapidity dependent average transverse 
mass (at midrapidity, % = my, see ref. [ESj f° r details), and T x and T y are direction (x 
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and y) dependent slope parameters: 

Here, Xf and Yf gives the transverse expansion rate of the fireball at freeze-out, and a = 
(T —T s )/T s is its transverse temperature inhomogeneity, characterized by the temperature 
at its center T , and at its surface T s . The thermal and collective contributions can be 
made more transparent by replacing p T with the transverse rapidity i/t = 0.51og((mT + 
Pt)/(^t —Pt)), to give the following approximate scaling law: 



k\ 2 ( -i i ^Tq ks ( Tq 



2 



T \ ki m ki \mj J 

where &i,2,3 are mass (m) dependent constants. Equation [7| shows that v 2 should scale as 
~ k m x my\ for different particle species or flavor. Hereafter, we define the scaled variable 
z/jf = k m x my\. The higher harmonic V4 ~ \v\ + k^y?, can be similarly obtained from 
Eq. It is also easy to show that v 2 scales with eccentricity and should be independent of 
the mass (size) of the colliding system for the same eccentricity. In what follows, we test 
whether or not such scaling is evidenced by RHIC data. The extent to which such scaling 
holds, gives an indication of the applicability and validity of perfect fluid hydrodynamics. 

3.1. Eccentricity scaling & system size dependence 
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Figure 20. Doubly differential anisotropy Figure 21. Scaled doubly differential 
V2(centrality ) pT) vs. Pt- Preliminary anisotropy V2(centrality,pT)/v2(centrality) 
PHENIX data are shown for several central- vs. pt- Preliminary PHENIX data are shown 
ity selections for Au+Au collisions [HZj. for several centrality selections for Au+Au 

and Cu+Cu collisions [ET7j. 

Fig. QUI show the differential V2(centrality,pT)- ) for mesons for several centralities. It 
shows the very well known dependence of v 2 on p T and centrality. Fig. I2T1 show v 2 results 
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for Cu+Cu and Au+Au scaled by the ^-integrated v 2 value obtained for each centrality 
selection. We note here that the ^-integrated flow is monotonic and linearly related 
to the eccentricity over a broad range of centralities. The scaled data for Au+Au and 
Cu+Cu shown in Fig. EJ are clearly compatible with the predicted eccentricity scaling. 
The fact that Cu+Cu and Au+Au scale to the same value indicates that the scaled v 2 
values are independent of system size as predicted. 

3.2. Particle flavor, pseudo-rapidity & higher harmonic scaling 
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Figure 22. Differential anisotropy v 2 (pt) for Figure 23. Differential anisotropy v 2 {y^ T s ) fo r 
protons, pions and kaons vs. pt [EZj- protons, pions and kaons vs. y^ s \ see text. 

Fig. E2l shows a comparison of the differential anisotropy v 2 (pt) for protons, kaons and 
pions obtained for the centrality selection 5-30%. The well known and rather characteristic 
flavor dependence of v 2 is clearly exhibited in the figure. If this aspect of the fine structure 
of v 2 has a hydrodynamic origin then the prediction is that it should scale with the variable 
y^ s = k m x my\. Fig. shows v 2 data for the same particle species scaled by the fine 
structure scaling variable y^ s . The results indicate scaling over a relatively broad range 
in yjf. It is interesting to note here that scaling appears to break in the y^ range where 
quark number scaling seem to work On the other hand, the latter is based on the 

coalescence of quarks which are flowing. Thus, hydrodynamic flow appears to dominate 
the dynamics for both low and intermediate pt particles. 

Fig. YIM shows v 2 data for a much larger selection of particle species (from both PHENIX 
and STAR) scaled by the fine structure scaling variable y^ s . Here again, the results indi- 
cate scaling over a relatively broad range in y^ s . The observed flavor scaling in these data 
provides rather strong evidence that the observed anisotropy is derived from a hydrody- 
namic origin. This conclusion is also supported by the good agreement obtained between 
the measured V4 and the scaled values {y\ ~ \v\ + k m y^) shown in Fig. ESI 

An excitation function of the pseudo-rapidity dependence of v 2 also serves as an impor- 
tant test for the predictions of perfect fluid hydrodynamics. Fig. EEI shows results from a 
recent analysis of PHOBOS data [EZj- The theoretical curves shown in the figure, give 
strong indications that these data are consistent with the analytic predictions of perfect 
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The universal scaling discussed above, clearly speaks to the validity of perfect fluid 
hydrodynamics over a broad range of the observed data. However, such scaling does not 
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provide explicit evidence for the degrees of freedom in the flowing matter. For intermediate 
values of pt ) hydrodynamic scaling breaks down (cf. Fig. EBJ) and the degrees of freedom 
can reveal themselves. Fig. E3 shows the results from a recent test for quark number 
scaling of v 2 by the STAR collaboration [ 07j . The lower panel of the figure illustrates 
the deviations from the quark coalescence ansatz used to scale the data shown in the top 
panel. The deviations are largest at low where hydrodynamic scaling was shown to 
work best. For intermediate p^s (where the fluid dynamical picture breaks down) the 
deviations are quite small, suggesting that the active and relevant degrees of freedom are 
those of constituent or valence quarks. 

4. Conical flow revisited 

Objects moving at supersonic speeds create "conical flow" behind the shock front they 
create. Such flow were initially conjectured to occur in cold nuclear matter [EBj. However, 
cold nuclear matter was later found to be too dilute and dissipative to sustain conical flow. 
The recent discovery of jet quenching at RHIC [EHj suggests that jets created in heavy 
ion collisions deposit a large fraction, if not all, of their energy into the produced matter. 
If this matter is indeed a strongly coupled Quark-Gluon Plasma (sQGP) having a small 
viscosity, then the current conjecture is that this energy should propagate outward in the 
form of conical flow as illustrated in Fig. El [ 9j. Extensive searches for such flow, via two- 
and three-particle correlation functions, are currently underway [ IfiHl I5T] . 




Figure 28. Hadron-hadron- Figure 29. Hadron-meson- Figure 30. Hadron-baryon- 
hadron correlation function [ meson correlation function [ baryon correlation function [ 

EOj. EDI EOj. 

Figs. EHI - EDI show three-particle Ac/) jet correlation surfaces (A(j)i^ vs. A^i^) for 
a trigger hadron from the range 2.5 < Pt <4.0 GeV/c (hadron #1) and two associated 
hadrons from the range 1.0 < Pt <2.5 GeV/c (hadron #2 and #3) [ 61J. Results are shown 
for the centrality selection 10-20 %. They show a strong dependence on the flavor (PID) 
of the associated particle and clearly do not follow the expected patterns for a "normal 
jet" [ 61j. Instead, they provide rather compelling evidence for strong modification of 
the away-side jet. Further detailed quantitative investigations are required to firm up 
whether or not the tantalizing qualitative features evidenced in Figs. EHI and EH confirm 
the presence of conical flow. 
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5. Epilogue 

In summary, the gross elliptic flow patterns observed at RHIC provide compelling ev- 
idence for the production and rapid thermalization (less than 1 fm/c after impact) of 
nuclear matter having an energy density well in excess of an order of magnitude above 
the critical value required for deconfinement. The universal and fine structure scaling 
properties of this flow, indicate that this new form of matter behaves like a perfect liquid. 
The quark number scaling observed at intermediate pts (when the fluid dynamical picture 
breaks down), indicates that the relevant degrees of freedom are those of constituent or 
valence quarks. These properties are consistent with those of a strongly coupled plasma 
with essentially perfect fluid-like properties. Systematic experimental and theoretical ef- 
forts are currently underway to quantitatively pin down the thermalization time, EOS 
and transport properties of this plasma. 
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